Q-dependence of the inelastic neutron scattering cross section for molecular spin 

clusters with high molecular symmetry 
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For powder samples of polynuclear metal complexes the dependence of the inelastic neutron 
scattering intensity on the momentum transfer Q is known to be described by a combination of 
so called interference terms. They reflect the interplay between the geometrical structure of the 
compound and the spatial properties of the wave functions involved in the transition. In this work, 
it is shown that the Q-dependence is strongly interrelated with the molecular symmetry of molecular 
nanomagnets, and, if the molecular symmetry is high enough, is actually completely determined by 
it. A general formalism connecting spatial symmetry and interference terms is developed. The 
arguments are detailed for cyclic spin clusters, as experimentally realized by e.g. the octanuclear 
molecular wheel Crg, and the star like tetranuclear cluster Fe4. 



I. INTRODUCTION 

Molecular nanomagnets, inorganic complexes with 
tens of magnetic metal centers linked by organic ligands, 
are a new class of magnetic materials exhibiting fascinat- 
ing quantum properties. For instance, quantum tunnel- 
ing of the magnetization has been observed in clusters 
like Mni2 and Fesji now called single molecule magnets. 
Another class of molecular nanomagnets is established 
by the molecular wheels as represented by the prototype 
Feio^ These antiferromagnetic cyclic clusters show pro- 
nounced steps in the magnetization curve at low temper- 
ature signaling quantum size effects^^ Furthermore, the 
possibility of coherent tunneling of the Neel-vector has 
been suggested in these compounds. ■ Yet other systems 
like molecular gridsr&LS tetranuclear star like clusters^ 
the cluster V 15 ^ or the magnetic Keplerate Mo72Fe3(jii 
attracted much interest due to their peculiar magnetic 
properties. 

The magnetism of these polynuclear complexes is, in 
principle, described by a microscopic spin Hamiltonian 
which to a very good first approximation consists of terms 
representing isotropic Heisenberg exchange interactions 
and dipole-dipole interactions among different spin cen- 
ters within a molecule, and terms due to the ligand field 
and Zeeman interaction for each single ioniZ^ Inter- 
molecular interactions are usually very small and can be 
neglected. The isotropic Heisenberg terms are dominant 
in general and these clusters may be thus described as 
(antiferromagnetic) Heisenberg spin clusters with weak 
magnetic anisotropy. 

For a detailed understanding of the properties of 
molecular nanomagnets it is of great importance to be 
able to study both the effects due to the Heisenberg inter- 
action and anisotropic terms. Although rarely achieved 
in praxis for these rather large molecules, the ultimate 
goal would be to determine all relevant parameters of the 
microscopic spin Hamiltonian like exchange constants, 
single ion zero-field-splittings, etc.. Besides the obvi- 
ous relevance in understanding the physics, the results 
are of great help also for synthetic chemists as they al- 



low to establish magneto-structural correlations which in 
turn provide ideas for a controlled improvement of these 
systems S 

Inelastic neutron scattering (INS) has been demon- 
strated to be a powerful experimental tool in this 
regardii^i^iSiil Being a spectroscopic technique, INS 
provides direct access to energy splittings. The INS 
cross section is controlled by the favorable selection rule 
AS = 0, ±1. INS thus allows to detect the splitting of 
individual spin multiplets, as EPR (AS = 0), but also 
the splitting produced by the magnetic interactions. As 
a result, INS provides a very straightforward determina- 
tion of exchange parameters. It requires, however, large 
amounts of sample and, frequently, deuteration. Accord- 
ingly, only powder samples of molecular nanomagnets 
were investigated so far. 

As a unique feature, INS additionally allows to mea- 
sure the dependence of the INS intensity on the scattering 
vector Q. In this work only the situation of powder sam- 
ples will be considered, i.e. the INS cross section has to be 
averaged in Q space and the scattering intensity will de- 
pend only on the momentum transfer ( t ) JJJi On the one 
hand, the Q-dependence allows to distinguish magnetic 
excitations from vibrational excitations^ On the other 
hand, and more importantly, the Q-dependence allows 
discriminating clearly between various types of magnetic 
transitions Thus, the observed Q-dependence of 
a transition can be used in addition to its energy for a 
spectroscopic assignment removing ambiguities in many 
cases. This advantage has been explored in several works, 
a particular illustrative example is given by the tetranu- 
clear cluster [Co 4 (H 2 0)2(PW9 34)2] for which the Q- 
dependence was key to establish the correct models 

As for the origin of the Q-dependence, it is well known 
that it depends sensitively on the wave functions of 
the states involved in a transitions^ In the formula 
for the INS cross section the matrix elements for in- 
dividual spin centers are correlated to their geometri- 
cal arrangement and the momentum transfer via a term 
exp [«Q(Rj — Rj)]. Here, denotes the position vector 
of the z-th metal center. This gives rise to characteris- 
tic interference terms in the INS cross section for powder 
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samples^ This observation was expressed in Ref. 0] 
loosely as 



II. INS CROSS SECTION FOR POWDER 
SAMPLES 



The differential neutron scattering cross section i sH4£ 



I(Q) oc F 2 (Q) [combination of interference terms] , (1) 



where I(Q) denotes the averaged INS intensity and 
F(Q) the magnetic form factor. It is thus clear, and well 
understood, that the Q-dependence basically senses the 
spatial properties of the involved wave functions. How- 
ever, to the best of the authors knowledge, beyond this 
basic understanding the significance of Eq. remained 
unexplored for larger spin systems like molecular nano- 
magnets. 

For a limited number of small sys- 
tems the Q-dependence could be calculated 
analytically, 14 i 19 i 21 i 22 i 23 i 24 thus detailing Eq. Q for 
these cases. For larger, or more complex systems one 
has to resort to numerical procedures^SiSSiSSiS! and a 
general program has been developed^ However, the 
"black box" character of this approach did not yield 
insight into the relation for large spin systems between 
geometrical structure and wave functions as expressed 
by the interference terms. 

This work aims at showing how the Q-dependence is 
influenced by the molecular symmetry of the spin clus- 
ter and in particular, if the molecular symmetry is high 
enough, that it is actually predetermined by it. The ar- 
guments will be detailed for the particular cases of the 
cyclic spin clusters with focus on the molecular wheel 
Cr§rS and the tetranuclear star like cluster Fe^ The 
molecular wheel Crs and star like cluster Fe4 were chosen 
as examples as they both exhibit a very high molecular 
symmetry and detailed INS experiments were performed 
recently^S^i*^ Additionally, the Heisenberg interaction 
is much stronger than the anisotropic terms in both com- 
pounds. INS measurements on powder samples thus pro- 
vide an excellent view on their internal spin structure. 
Several transitions between different spin multiplets were 
observed experimentally and only these are of interest 
here. As a main result it will be demonstrated that for 
such highly symmetric systems the Q~dependence actu- 
ally can be used to determine experimentally the spatial 
quantum numbers of the states involved. The molecular 
wheel Crs is a particularly illustrative example in this 
respect. 

In the next section, first the general equation for the 
INS cross section of powder samples will be derived, cor- 
recting an earlier resulti 2 ^ Then, in section HTTl the general 
framework based on standard group theoretical proce- 
dures will be developed. In sections llVl and Ivl the molec- 
ular wheel Crg and the star like cluster Fe4 will be dis- 
cussed. The work concludes with section fvTI 
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where C(Q,T) = (^e 2 /m e c 2 )(k' /k) exp[—2W(Q, T)], 
(3 = l/(fcsT), Z(T) the partition function, and 



J„ m (Q) = ^^(Q)F i (Q)e i Q' R - 
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)(n\S ia \m)(m\S j i 3 \n). 



(3) 



In this equation, Fi(Q) is the magnetic form factor 
of the i-th spin center, a, [3 — x,y, z, Q = k' — k is the 
transferred momentum, and Ry = R; — Rj with the spin 
position vectors R^. In the following, the abbreviations 
la = Qa/Q and Si a for the matrix elements of Si a are 
used. In an ordered product Si a Sjf3, the first term always 
has to be understood as a matrix element between (n| and 
|m) and the second as between (m| and |n). 

The cross section for a powder sample is obtained by 
averaging Eq. (J2J over all directions of Q, i.e. one has to 
calculate the expression 
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R * J 5^(^ Q /3 — lalp)SiaSjp. (4) 
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This is conveniently done by resorting to the calculus of 
irreducible spherical tensors^ The sum over a and (3 can 
be rearranged to J2 a ^ia^ja - (Y, a l*Siod{Y!,plpSjp)- 
The first term in this expression is easily averaged since 
it does not depend on the orientation of Q, yielding 
jo(QRij)Si ■ Sj. Here, jk is the spherical Bessel function 
of order k. The second term is calculated by express- 
ing the scalar product ~}2 a l a Sia by spherical tensors^ 
according to 



^^a = E T i 1) *(Q) T i 1) (^)- 



(5) 



Tq k ^ (O) denotes the q-th. component of the irreducible 
spherical tensor of degree k with respect to the vector O. 
They are proportional to the spherical harmonics Y^(0). 
Expanding the exponential in spherical harmonics and 
using the coupling rule for spherical tensors one obtains 
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Here, (...) denotes a Wigner-3j symbol. The term 
for L = gives \jo{QR l j)S l ■ Sj. For I nm {Q) = 
J dfi/47rJ„ m (Q) one finally obtains 



WQ) =J2F*(Q)F j (Q)(~j (QR ij )S i ■ 8j 
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The L = 2 term could be put also into the more com- 
pact form j 2 (Qi?y )T( 2 )(Ry ) • T( 2 ) (SiS-,) using Eq. J5J. 

Several special cases shall be discussed. It is interesting 
to consider the limit Q — > 0. Since jo(a; — + 0) » 1 + C(a; 2 ) 
and j2(a; — ► 0) w x 2 /5, one simply obtains 

4m(Q -> 0) = ~Y / F?(Q)F j (Q)S i ■ Sj + 0(Q 2 ). (8) 

Next, for an isotropic spin cluster the L = 2 term in 
Eq. (JJJ obviously vanishes. If one further numbers the 
different possible values of R+j by p and writes R p , the 
result assumes the convenient form 
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(Q,T) = J2^ m iQ,T)j (QRp)- 



(9) 



Accordingly, the scattering intensity is sort of an ex- 
pansion in jo(QRp). In principle, since the R p are known, 
the coefficients I p can be determined by fitting the (in- 
dependence of the scattering intensity to Eq. ©. Equa- 
tion iJS} precises Eq. for the case of isotropic systems. 

Finally, one may consider the case where only the 
components with a = (3 contribute. This holds if 
(Si t0l Sjp + SipSja) = for all a ^ j3. However, it is 
clear that this condition is not fulfilled in general. One 
necessary condition is that the quantization axis for the 
spin operators coincide with a magnetic main axis. For 
uniaxial and isotropic systems this is sufficient to ensure 
the above condition. However, for a biaxial magnetic 
system this is not sufficient in general since in the neu- 
tron scattering cross section the local terms Si a Sjp enter 
while " biaxial" is a property of the whole spin cluster as 
an entity. Thus, for spin systems with magnetic symme- 
try lower than uniaxial the validity of this simplification 
has to be checked carefully. For an uniaxial system only 
the M — term contributes, resulting in 



I^(Q) = ^F*(Q)Fj(Qy, 



h(QRij) - -j 2 (Q%)C 2 (R y ) 



(s ix s jx 



SiySjy 



+ [jo(QR ij )+32(QRij)Ci(R ij )] SiJjzllO) 
with Cg(Ry) = [ZiRi^/Riof - 1] /2. 



III. INS CROSS SECTION AND SPATIAL 
SYMMETRY 

The spatial symmetry of the molecule leads to an in- 
variance of the spin Hamiltonian upon certain permuta- 
tions of the spin centersi 34 i 35 The group of all these per- 
mutations is denoted as Q 1 and an element of this group 
as P. It should be noted that Q and the point group of 
the molecule are not equivalent in general. This symme- 
try of the spin Hamiltonian is accordingly denoted here as 
spin permutational symmetry^ Its effects are exploited 
by applying the standard results of group theoryi 33 i 34 

The eigenstates |n) of the spin Hamiltonian can be 
classified by the irreducible representations (IRs) of Q 
and are written now as \rk\i). They transform accord- 
ing to 0{P)\Tkp.) =E A r ^(- p )k&/i). 0{P) denotes the 
operator associated with the permutation P. Further- 
more, 0(P)Si a O(P~ 1 ) = Tji(P)Sj a - Here, no sum over 
j arises on the r.h.s since P is a permutation, i.e. Pi Pj 
for i j. Accordingly, each row or column, respectively, 
of the matrix T(P) has only exactly one nonzero entry 
which is equal to 1. With these relationships one confirms 



E 



{r'k'^S^rkp) = 



(11) 



where Pu is a symmetry element which maps the j- 



th spin exactly onto the i-th spin, i.e. Pjij 



Equa- 



tion (|llfl describes the restrictions imposed by the molec- 
ular symmetry on the interference effects in inelastic neu- 
tron scattering, and is thus central to this work. 

The spin centers can be divided into classes such that 
all the spin centers of one class are related by the permu- 
tations of the group Q. More precisely, a class is a set of 
spin centers which transform one into another under the 
operations of Q . The different classes will be numbered 
by 7. By construction, each spin center is member of 
exactly one class. Therefore, 



E = EE 



(12) 



7 i€7 



From each class, one spin center is chosen (arbitrar- 
ily) and called the pivotal center of this class. It will 
be also indicated by 7. Equation (|11|) then essentially 
provides a relation between the matrix element of an ar- 
bitrary spin center with the matrix element of its piv- 
otal center. This imposes a certain structure on the ex- 
pression F *(Q) F j(Q)f(Q,^4j)Si a Sjp, which, with 
/(Q,Rjj) appropriately chosen, is the part of the cross 
section relevant here. This idea can be used to work 
out general expressions based on the generalized Wigner- 
Eckart theoremi 3 ^ 

For spin permutational groups with only one- 
dimensional IRs, writing the eigen states as \rk), Eq. 
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simplifies to { T 'k'\S jP \Tk) = (T , k'\S t0 \Tk)x k ' k (P jl )- 
Whereby, \ kk ' (P) = X k *(P)x k '(P) and X k being the 
character of the fc-th 1R. Together with Eq. I|12() and af- 
ter some rearrangement of terms one obtains the central 
result of this work: 

ij 

Y / S ia S^F;(Q)F^Q)f kk '(Q) (13a) 

77 

f kk '(Q) = ^^/(Q,R y )x fcfc '(P 47 )x fc ' fc (^)- (13b) 

The Q-dependence of the neutron scattering cross sec- 
tion is given by the interference terms f kk (Q) which can 
be calculated without knowledge of the wave functions, 
i.e. are completely governed by the spatial symmetry 
properties. These equations, being a precise statement of 
Eq. 0}, become the more useful the smaller the number 
of classes 7 produced by the spin permutational symme- 
try is. 

IV. CYCLIC SPIN CLUSTERS 

The full spin permutational symmetry group of a cyclic 
spin cluster is I? at, where N is the number of spin cen- 
ters. However, the spin levels may be also classified by 
the IRs of the subgroup Cm reflecting the translational 
invariance^ This effectively introduces the shift quan- 
tum number q = 0, . . . , N — 1 via the shift operator T: 
T\rq) = e~ tq2n ' N \Tq). The group elements may be writ- 
ten as T n with n = 0, . . . , N — 1 and the characters are 
X (?)(T") = e -^ni^lN _ The fuU symmetry group Dn has 

one-dimensional and two-dimensional IRs. Spin levels 
with q = and q — N/2 belong to one-dimensional IRs 
of Dn, while for q 7^ 0,N/2 the two states with q and 
N — q belong to a two-dimensional IR. The latter are thus 
degenerate. 

The interference terms f kk (Q) can now be calculated. 
For each spin center there is always one symmetry el- 
ement T n which connects it to any other spin center. 
Thus, there is only one class 7: 7 = {1, . . . , N}. Spin 
center 1 is chosen as pivotal center. For P^ one es- 
tablishes that Pi 7 = Pn = T 1 ^ 1 , and accordingly that 

X qq '(Pij)x q ' q (Pj-y) = e^ q - q '^- l ^/ N . Introducing the 
" distance" n by j = i + n yields 

N-l N 

f qq '(Q) = ^J]/(q I H i ,H»)^ (M ' ) ". (14) 

n=0 i=l 

The cyclic spin clusters synthesized so far exhibit at 
best a molecular Sn symmetry axis, as determined by 
x-ray crystallographyi^^LSi The magnetic anisotropy 
of these systems is thus expected to be strictly uniax- 
ial, as is supported by experiments. 39 '— The position 



vectors point to above and below the plane of the 
molecule perpendicular to the main symmetry axis. How- 
ever, even with this effect /(Q, His+ n ) becomes inde- 
pendent on i for powder samples, i.e. Ri^+ n — > R n and 
C^(R M+ „) -> C$(Rn) [see Eq. lfTU)l]. Thus, for powder 
samples of these highly symmetric, cyclic spin clusters 
the interference terms are given by 

N-l 

f qq '(Q) =NJ2 /(Q,Rn)e^ (9 - 9 ' )n . (15) 

n=0 

The averaged INS cross section for the transition 
\rq) — > \r'q') is then simply proportional to 

^rv(Q) = \nF 2 {Q){S 1 • S 1 /« 9 ,(Q) + 
i (2S lz S lz - S lx Sxx - SiySiy) f m , (Q)} (16a) 
with 

n=0 
N-l 

flAQ) = E MQR n )c 2 (R n y^-i>. (16c) 

Apart from two normalizing factors, the Q-dependence 
of the INS intensity of uniaxial cyclic clusters is analyt- 
ically determined by their spatial symmetry properties. 
The interference terms f qq , (Q) and f 2 , (Q) do not de- 
pend on t and r'; they are completely determined by the 
spatial quantum numbers. 

The following general result is noteworthy. In the limit 
Q -> 0, Eq. leads to 

AT-l 

I rqT , q ,(Q) = -NF 2 {Q)SlY,^ [q ~ q ' )n + 0{Q 2 ) 

n=0 

oc 5 qq , +0{Q 2 ). (17) 

Therefore, comparing with Eq. JSJl, the spatial symme- 
try of the cyclic spin cluster enforces the INS intensity 
to become strictly zero for Q = if q ^ q\ i.e. if the 
spatial quantum numbers of the states involved are dif- 
ferent. In this case, the intensity actually approaches 
zero as I(Q) oc Q 2 . Equation (|17fl might have important 
experimental implications. Recently it became an inten- 
sively discussed question whether cyclic clusters like the 
ferric wheels are accurately described by a spin Hamilto- 
nian with perfect cyclic symmetry, or whether additional 
terms with less symmetry are of importance4i^£ In prin- 
ciple, Eq. I|17[) provides a powerful experimental approach 
to unravel this question since if symmetry breaking terms 
are relevant the INS intensity would not drop to zero with 
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V2Ro, and 



Eq. I|16(l shall be worked out explicitly for the case 
of an octanuclear ring with predominantly antiferromag- 
netic next neighbor interactions, as appropriate for the so 
called Crg cluster. 29 This is motivated by recent, very de- 
tailed and accurate INS measurements on this system^ 
The Q-dependence of the INS intensity has been mea- 
sured earlier also on a hexanuclear iron(III) cluster, but 
has been found afterwards to be less reliable due to in- 
strumental artifactsi 3 ^ 3 . 

For an octanuclear ring with molecular S4 sym- 
metry, the position vectors are simply Rj = 
Ro (cos(«7r/4), sin(i7r/4), (— l) l R z ), where Ro is the "ra- 
dius" of the ring (Ro = 4.427 A for Cr 8 ) and R z deter- 
mines by how much the ions lie above or below the plane 
of the ring. Since R z <§C Ro, it shall be neglected here 
(but can be incorporated straightforwardly if felt to be 

required). For R n then holds R\ — Rj — 

R2 = Rq = V / 2i?o, R3 = R5 = V? 
i? 4 = 2Ro. Also, it follows that C^(R„) = -1/2. This is 
actually valid for the weaker condition R 2 Z <C Rq. 

For finite rings with antiferromagnetic Heisenberg in- 
teractions it is well established that, as function of the 
total spin quantum number S, the lowest lying states 
form a set of rotational bands,— The energies of the 
states of a particular rotational band follow the Lande 
rule E(S) cx S(S + 1) as for a rigid rotatoria The 
set of bands can be distinguished into L- and L-band: 
The L-band consist of the states for which q toggles as 
function of S between q = and q = N/2, the L-band 
embraces the lowest lying states with q ^ 0, N/2, which 
are essentially spin waves in character— 4& The energies 
of the states of these bands can be approximated excel- 
lently by E(S, q) = ±AS(S + 1) + e(q). A characterizes 
the "curvature" of the bands and coincides with the en- 
ergy gap between ground and first excited state. e(q) 
measures the "offset" of a particular band, as classified 
by q. Whereby, e(q) ~ const] sm(2irq/N) Impor- 
tantly, the L-band is well separated in energy from the 
higher lying L-bands, i.e. A <g; e(q) for q 7^ 0, N/2. 

In actual INS experiments only the transitions between 
states within the L-band and those from the L-band into 
the L-band are relevant since A is typically on the order 
of meV. For transitions within the L-band q — q 1 = N/2 
holds; and for transitions from the L-band into the E- 
band q — q 1 7^ 0, N/2 (values for q should always be un- 
derstood modulo N). 

First, the L = contribution in Eq. fTfif. i.e. fq q i{Q) 
will be considered. For Crg, q — q' = 4 for transitions 
between states of the L-band, directly yielding 



a) 



J 0,. 



l + 2j2(-l) n j (QR n )+j (QR 4 ) (18a) 



for the interference term. For transitions from a state 
of the L-band to states of the L-band one additionally 
has to take into account the degeneracy of the states with 
q and N — q due to the overall Dn symmetry. However, 
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FIG. 1: (a) Q-dependence of the interference terms /q 4 , /q x , 
/o 2 an d /o,3 a s given in Eq. 118H for an octanuclear cyclic ring, 
(b) Dependence of the inelastic neutron scattering intensity 
as function of momentum transfer Q for the three observable 
transitions at low temperatures of an octanuclear antiferro- 
magnetic Heisenberg ring with spin-3/2 (see text). For both 
panels Ro = 4.427 A was used as appropriate for the cyclic 
molecular cluster Cr8. 



inspection of Eq. I|15fl shows that the interference terms 
are equivalent for q' and N — q' (and of course analogous 
for q). One finds 



f° - 

Jo.i — 

f° - 

J 0,2 — 

f° - 

Jo, 3 — 



1 + V2 \jo(QRi) - MQR 3 )} - j (Qi?4)(i8b) 
l-2j (QR 2 )+h(QRi) (18c) 
1 - V2 \j (QRi) - MQR3)] - MQRiXl&d) 



These functions are plotted in Fig. 1(a) with Rq cho- 
sen as appropriate for the Crs cluster. To compare with 
experiment one additionally has to consider the quasi 
degeneracy of the spin levels with q — 1,7 and q = 3, 5 
because of the | sin(27rg/A)| like dependence of e(q) on 
q. Since this quasi degeneracy cannot be resolved exper- 
imentally, one should sum these contributions for a com- 
parison with observable peaks. This requires a knowl- 
edge of the oscillator strengths = ~}2, a {Tq\Si a \T' q') 2 



G 



TABLE I: §1 = T, a {Tq\S la \r'q'} 2 for the possible transitions 
from the ground state (S = 0, q = 0) to the L-band (S = 
1, q = 4) and to the S-band (S = 1, q' = 1, 2, 3, 5, 6, 7) for an 
octanuclear antiferromagnetic spin-3/2 Heisenberg ring. 



q'= 


1,7 


2,6 


3,5 


4 




0.06802 


0.1870 


0.4359 


2.316 



for each of the four contributing transitions. Due to sym- 
metry reasons the oscillator strengths for q' and N — q' 
are equivalent. The dependence on q is approximated 
by y/1 - cos(q2ir/N)/l + cos(g27r /N)£L but for precise 
results the oscillator strengths should be calculated, e.g. 
numerically. For the possible transitions starting at the 
S = ground state, they are given in Table [I] 

Combining Eqs. I|16l) and (|18|) and using Table |U the 
Q-dependence of the INS intensity for the three transi- 
tions as observable for Crs at low temperatures can be 
expressed analytically and is shown in Fig. 1(b) in units 
of 2/3NF 2 (Q). These curves were found to be in excel- 
lent agreement with experiment^ 

Since Cg(R„) = —1/2, the calculation of fg q /{Q) pro- 
ceeds in exactly the same way as for /° , (Q). One obtains 
basically Eqs. I|18l) multiplied by —1/2 and with jo(QR n ) 
replaced by j2(QRn)- Thus, the magnetic anisotropy in 
the cyclic spin clusters mixes in some j% contributions, 
but since the anisotropy in general is weak compared to 
the antiferromagnetic coupling, the Q-dependence of the 
INS intensity is rather well described by f qq , (Q) alone. 



V. TETRANUCLEAR STARLIKE SPIN 
CLUSTERS 

For further illustration of the use of Eq. (|13|) . the 
tetranuclear starlike clusters will be considered^ In these 
clusters a central spin is surrounded by three spins which 
form an almost perfect equilateral triangle. Recently, 
such systems incorporating four iron(III) ions were inves- 
tigated and attracted considerable interest as they were 
found to behave as single molecule magnets and show the 
still rare phenomenon of resonant quantum tunneling. 9 
For one member of this class of compounds the INS spec- 
tra has been measured recentlyi^iLi^i The Q-dependence 
has not yet been measured reliably^ but can be expected 
to be investigated with better accuracy in the near future. 

The actual symmetry of these molecules is less than 
threefoldASi They exhibit a twofold symmetry axis 
which passes through the central spin and one of the sur- 
rounding spins. The magnetic anisotropy, accordingly, 
is biaxial, but much smaller than the antiferromagnetic 
coupling between central spin and surrounding spins. It 
shall thus be neglected here. This essentially corresponds 
to neglecting the L = 2 terms in Eq. J7J, similar to the 
above discussion for the cyclic clusters. Also, the geomet- 
rical deviation from threefold symmetry is very small. It 



is thus a very good approximation to take the spin per- 
mutational symmetry as D%. 

The arguments to calculate the interference terms are 
very similar to those for the cyclic spin clusters. The 
new feature which arises here is that the spins divide 
into two classes 7. Again, the subgroup C3 will be con- 
sidered. Since there is no group element T n (n = 0, 1, 2) 
which maps the central spin (numbered by 4 here) onto 
one of the surrounding spins, one obtains the two classes 
{1,2,3} and {4}. It should be noted here, that one 
obtains exactly the same classes even for the full sym- 
metry group -D3. The Q-dependence is determined by 
the interference terms ffl (Q), /f| (Q), fH (Q), and 
fit (Q) m Eq. I|13a|) . In the present case, however, 

f?l(Q) = fH'(Q) and fii(Q) = F 2 (Q)Si The av- 
eraged INS cross section is then 



Si • SxfliiQ) + Si • S 4 f?!(Q) + s 4 • s 



(19a) 



with 



fn (Q) = 3 + 6MQR12) cos ( —(q - q') ) (19b) 



2tt 



ffi (Q) = 3j (QRu)6 qq , 



(19c) 



Similar to the cyclic clusters, spin levels with q = 
belong to an one-dimensional IR and that with q = 1, 2 
to the two-dimensional IR of D 3 . The latter are thus 
degenerate and one should add the respective INS inten- 
sities. 

In the starlike clusters, the antiferromagnetic coupling 
between the spin-5/2 iron(III) ions leads to a S = 5 
ground stat o 9 i 30 belonging to q = 0. The first excited 
states, approximately 80 K above the ground state, are 
made up of two S — 4 levels 9 *^ belonging to q = 1,2. 
This excitation has been observed in experiment i^i Its 
INS intensity is calculated as 



Io,i2(Q) = ^F 2 (Qy, 



3{2Si-S 1 [l-j (Qi? 12 )] + S4-S 4 } 



(20) 



since, as for the cyclic clusters, the interference terms 
remain unchanged by q' — > — q' . Up to a constant, the 
Q-dependence is here again determined by the spatial 
symmetry properties of the molecule. 

Equation (|20[1 is expected to describe the Q- 
dependence of the observed S' = 5^2xS' = 4 transition 
well, not only because of the reasons given above, but 
also because the only structural element which enters is 
R\2- This implies a strategy to use the Q-dependence as 
a check for effects which break the D3 spin permutational 
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symmetry (the Q — > argument is not working here be- 
cause of the Q independent contribution fH ). Such ef- 
fects lead in particular to a contribution proportional to 
io(Q-Ru) which would be absent otherwise [Eq. J20J)]- Of 
course, L = 2 contributions also will differ, but they re- 
main small as the magnetic ansiotropy is much smaller 
than the coupling. 

Transitions from the ground state to even higher lying 
states or within the excited states should be very diffi- 
cult to observe experimentally and thus are not discussed 
here. But transitions within the S — 5 ground state mul- 
tiplet, corresponding to q = — > q' — 0, were observed^ 
The calculation of their Q-dependencies, however, is es- 
sentially just a repetition of the above considerations. 

VI. CONCLUSIONS 

In this work, it has been shown for the two cases of 
the cyclic and the star like clusters how the spatial sym- 
metry of the molecule predetermines interference terms 
and thereby the Q-dependence of the INS intensity. The 
symmetry elements of the spin permutational symmetry 
group Q impose relations among the matrix elements at 
different spin sites. The number of classes 7 allowed by 
the spin permutational symmetry turned out to be es- 
sential: The smaller the number of classes is, the fewer 
spin matrix elements remain uncorrelated, and the more 
effective the whole procedure becomes. The two cases 
discussed here are in a sense extreme as the number of 
classes 7 produced by the symmetry operations was very 
small, being even one for the cyclic spin cluster. 

The procedure of how to calculate the effects of the 
spatial symmetry has been worked out here only for the 
case of one-dimensional IRs. This might appear as a 
significant restriction, but it is sufficient for many prac- 
tical cases. For instance, in the cases presented here, the 
full symmetry group actually had two-dimensional IRs. 
This has been handled by first considering a subgroup 
with only one-dimensional IRs, calculating the interfer- 
ence terms for this situation, and then, in a second step, 
to take into account the degeneracies imposed by the full 
symmetry by simply adding the corresponding intensi- 
ties. In the present cases this was quite trivial as the 
interference terms were identical for the involved degen- 
erate states. In the general case, one has to work out 
the relationship between the involved matrix elements. 
However, since it is entirely determined by transforma- 
tion properties it can be calculated using the generalized 
Wigner-Eckhard theorem and the irreducible operators 
as appropriate for the symmetry group Qw^ 

This approach worked for the cyclic and star like clus- 
ters without loss of information about the Q-dependence 
because the symmetry elements of the subgroup were 
complete enough to produce the minimal number of 
classes. That is, exactly the same class structure as for 
the full symmetry group was obtained. Certainly, if this 
does not hold, one either has to accept additional un- 



determined matrix elements in the analytical expression 
for the Q-dependence, or has to start right away from 
the general, but somewhat involved, equation Eq. 

It is noteworthy that the outlined strategy works for 
quite a number of important cases. For instance, tetranu- 
clear clusters of the so called cubane structure are very 
frequent and constitutes an important class of tetranu- 
clear compounds. Actually, a significant number of 
tetranuclear clusters with a different topology (e.g. grids, 
squares, or chains) became available only recently due to 
advances in inorganic chemistry. The metal centers in 
these complexes approximate a tetrahedron, i.e. exhibit 
a rather high molecular symmetry. However, the C4 sub- 
group is sufficient to connect all spin centers by symmetry 
elements resulting in only one class 7 including all spin 
centers. Thus, these clusters can be handled in exactly 
the same way as the cyclic spin clusters. 

In many cases, the molecular cluster does not exhibit 
a high symmetry exactly, but approximates it with small 
distortions. The Bessel functions will be affected only 
slightly since deviations from the regular positions of the 
spin centers are smaller than the distances between the 
spin centers. Furthermore, the distorted structure also 
leads only to small changes in the parameters of the spin 
Hamiltonian, which typically have no effect on the en- 
ergy spectrum up to first order in perturbation theory^ 
As a result, the interference terms are actually rather in- 
sensitive to weak distortions from an assumed optimal 
symmetry. The Q-dependencies as calculated for the ap- 
proximate symmetry are then nonetheless useful, in par- 
ticular as the accuracy of experimental Q-dependencies 
is often limited. 

With regard to the behavior of the INS intensity for Q 
approaching zero, it has been found for the cyclic clusters 
that it drops to zero if states with different spatial spin 
quantum numbers are involved in the transition. Vice 
versa, it will approach a finite value only if the spatial 
spin quantum numbers coincide. This is a very general 
feature for clusters for which all spin centers are con- 
nected by elements of the symmetry group G, i.e. for 
which only one class 7 appears. With this reasoning in 
mind, the different behavior of the Q-dependencies of 
the transitions observed for e.g. the dimer Tb 2 Br^ _ be- 
comes immediately transparent^ - it is a result of the 
inversion symmetry of the dimer. These differences for 
Q — ► were used in this case with advantage to con- 
firm the assignment of observed peaks to theoretically 
expected transitions.— As is clear from the discussion of 
the tetranuclear star like cluster, this kind of argumen- 
tation does not apply to cases in which more than one 
class 7 arise. 
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